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Abstract 

This paper investigates the reconstruction of elastic Green's function from the cross-correlation 
of waves excited by random noise in the context of scattering theory. Using a general operator equa- 
tion, -the resolvent formula-, Green's function reconstruction is established when the noise sources 
satisfy an equipartition condition. In an inhomogeneous medium, the operator formalism leads to 
generalized forms of optical theorem involving the off-shell T-matrix of elastic waves, which describes 
scattering in the near-field. The role of temporal absorption in the formulation of the theorem is 
discussed. Previously established symmetry and reciprocity relations involving the on-shell T-matrix 
are recovered in the usual far-field and infinitesimal absorption limits. The theory is applied to a 
point scattering model for elastic waves. The T-matrix of the point scatterer incorporating all re- 
current scattering loops is obtained by a regularization procedure. The physical significance of the 
point scatterer is discussed. In particular this model satisfies the off-shell version of the generalized 
optical theorem. The link between equipartition and Green's function reconstruction in a scattering 
medium is discussed. 



1 Introduction 



Recently, there has been increasing interest in the use of coda and scattered waves to monitor temporal 
variations in dynamic structures such as volcanoes and faults. Such an approac h is known as coda 



wave interferometry in seismology and diffusive wave spectroscopy in acoustics (jSnieder et all 12002 



Snieder and PaeeL l2007h . While the early de velopments of the technique were based on the use of re- 
peating small earthquakes termed "doublets" ( Poupinet et ai , 1984, 2008|), the possibility to reconstruct 



Green's function f rom ambient seismic noise has given rise to a new technique termed passive image in- 
terfero metry (jSens-Schonfelder and Wegler . 200G; IWegler and Sens-Schoenfelderl 120071: IBrenguier et al 



l2008al lbh. In coda-wave interferometry, temporal changes in the medium are detected by comparing the 
coda portions of the cross-correlation function of ambient noise, that were computed at different times. 
Therefore, understanding the basic physical processes that allow the reconstruction of scattered waves 
from the cross-correlation of random noise sources is of interest in seismology and acoustics. 

The rec onstruction of G reen's function from field cross-correlations has a rich and vast history as 
reviewed by IShapirol ( 2011 ). Green's function retrieval in seismology and elastodynamics finds its roots 
in the fluctuation-dissipation theor em and in the theory of speckle correlations in optics. In recent years, 
the pioneering results described bv lShapiro ( 201lh have been extended to different system configurations: 
open or closed, and to different types of excitations: deterministic or random sources located on a surface 
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or distributed in a volume. For instance, general theorems on t he reconstruction of Green's function 
have been established for a variety of inhomogeneous, d issipative ( Wapenaar et al. . 20061 Snieder et al. . 
2007), and even non-reciprocal media (jWapenaarl I2006T ). An important ingredient of seismic interferom- 
etry is Green's theorem or it s elastodynamic generalization, i.e. the Rayleigh-Betti reciprocity theorem 
([Ramirez and Weglein , l2009h . The importance of Green's theorem and its gener alizations for the recon- 
struction of Green's function is well illustrated by I Wapenaar and Fokkem a (2006). Depending on the ex- 
citation mechanism -active or passive-, the acquisition geometry -surface or borehole measurements-, the 
nature of boundary conditions -radiation or reflection-, the type of Green's function employed -analytical, 
numerical, empirical-, a great variety of interferometric methods can be derived from the application of 
Green's theorem. A comprehensive review of applicati ons to exploration seismolog y, emphasizing the 
strengths and limitations of each method is provided by iRamirez and Weglein ( 20091) . 

More specifically, the connection between the cross-correlation of wavefields an d Green's functio n 
between two points in an inhomogeneous elastic half-space has been established by IWapenaarl (|2004l ). 
based on Betti's reciprocity th eorem. The reconstruction of Green 's function of a homogeneous elastic 
medium has been discussed by ISanchez-Sesma and Campillol (2006) from a different perspective. These 
authors find that a set of uncorrelated plane P and S waves verifying the equipartition relation, i.e. 
carrying the same amount of energy, allow the elastic Green's function to be reconstructed from the cross- 
correlation of nois e wavefields. This result wa s extended to the case of one cylindrical inclusion in a 2-D 
elastic medium bv lSanchez-Sesma et al. ( 20061) . These authors demonstrated the complete reconstruction 
of Green's function including scattered waves in the case of an illumination of the heterogeneous medium 
by an equipartition mixture of plane P and S waves. 

The case of scattering media has been studied by ISatol (|2009l 2010). He demonstrated the recon- 
struction of singly-scattered coda waves from random noise sources in an acoustic medium composed of 
point-like velocity perturbations. He considered both volumetric and surface distribution of noise sources. 
An important connection between energy conservation in scattering and the reconstruction of scattered 
waves in G reen's function obtained fr om cross-correlations of random noise sources has recently been put 
forward bv lSnieder et all (|2008ll2009[). Considering a scattering medium illuminated by noise sources dis- 
tributed on an enclosing surface JSnieder et al. ( 2008t ) showed that the reconstruction of Green's function 
of scalar waves is equivalent to a generalized optical theorem. Such theorems provide general symmetry 
relations that must be obeyed by the scatterin g amplitu d e, wh ich describes the amplitu de and phase 



relatio ns between incident and scattered waves. lLu e£ all (| 20 111 ) extended the method of lSnieder et all 
(2008) to the case of vector elastic waves, for a scat tering medium composed of a homogeneous and 
isotropic heterogeneity. ISnieder and Fleurv (l2010l ) and iMargerin and Satd (|201ll ) developed a theory for 
the reconstruction of scattered arrivals in Green's function obtained from the cross-correlation of noise 
signals in an acoustic medium, for surface and volumetric distributions of noise sources, respectively. 
These authors demonstrate the cancellation of spurious terms in the multiple scattering series expansion 
of Green's function cross-correlations by repeated application of the optical theorem. 

In this work, we analyze the case of a general -possibly anisotropic- heterogeneity embedded in an 
isotropic medium, illuminate d by a volumetric distr ibution of random, uncorr elated forces. In the same 
spirit as Snieder et all ( 20081 ). Lu et al. ( 2011 ). and Margerin and Satol ( 2011 ). we establish the equiva- 
lence between Green's function reconstruction and some symmetry relations pertaining to the scattering 
of elastic waves which we refer to as "generalized optical theorems" . Compared to our previous work 
on scalar waves, we first use a general operator formalism to establish in a concise way Green's function 
reconstruction theorem for a general -homogeneous or inhomogeneous- elastic medium. Next, we apply 
the theorem to Green's function of a scattering medium to establish general symmetry relations in the 
scattering process. An analysis based on scattering theory finally shows that the symmetry relations are 
extensions of generalized forms of the optical theorem for elastic waves. The central result of this paper 
-Equation (61)- contains as special cases all forms of optical theorems found in the literature. In view of 
future applications to a collection of scatterers we develop a simple point-like scattering model for elastic 
waves. The formulas ff80|) - ()81[) present an extension of the point-scattering model originally developed 
for scalar waves to the case of vector elastic waves. The physical interpretation of the model is discussed 
in details. 
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Figure 1: Problem setting. We consider the reconstruction of Green's function including the scattered 
waves from the cross-correlation of random wavefields recorded in a homogeneous isotropic elastic medium 
containing an arbitrary heterogeneity of typical dimension D (grey region). Receivers may be placed 
at arbitrary positions in the medium, including the near-field or the interior of the scatterer. The 
applied forces are stationary in time, randomly and uniformly oriented, and randomly and homogeneously 
distributed. The effective radius of the sources is R e = c s t with r the absorption time and c s the shear 
wavespeed. The effective radius is (much) larger than the size of the heterogeneity D and the typical 
distance between the receivers. 
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The problem setting is depicted in Figure [T] We consider an arbitrary heterogeneity with bounded 
support embedded in an elastic medium illuminated by a randomly homogeneous distribution of forces 
oriented in all possible space directions. By considering the correlation of the displacements recorded 
at two points A and B, we will show the possibility to reconstruct the elastic Green's tensor including 
the waves scattered by the heterogeneity. In a heterogeneous, slightly dissipative elastic medium, the 
wavefield generated by random forces f(x, t) satisfies the following equation: 

( Qf2 + ~ Qf- ) + M*) gf 2 d Xj C^ kl d XkUl (x,t) 

- 0, (i( ",,!.,<).,-. "/(x./i = fi(-K,t), (1) 

where po and C* fei = XodijSki + po($ik5ji +5u5jk) denote the density and elasticity tensor of the isotropic 
matrix. The summation over repeated indices is assumed. The second term on the left-hand side in- 
troduces a characteristic absorption time r. The deviations of the density and elastic properties from 
homogeneity are encapsulated in the terms dp = p — po and SCijki — C^i — C^ kl . No particular sym- 
metry properties, except the usual ones, are imposed on Cyjw. On the right-hand side of equation ([1]), 
we consider a randomly homogeneous and stationary distribution of forces with uniform random orienta- 
tions. In addition, we assume that the spatial correlation of these forces is much shorter than any other 
scale length of the problem, which implies that they are well described by a white noise process in space. 
The basic observable in seismology is the cross-correlation of signals averaged over an ensemble of noise 
sources: 
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T/2 



C ij (x B ,x A ;t)= lim -/ ((u i (x B ;t')uj(x>i;* / -t)*))dt', (2) 



T/2 



where the double brackets ((•}} denote the ensemble average. If the noise sources are statistically stationary, 
the basic quantity to be evaluated is the integrand on the right-hand side of equation ([2]). Note that the 
definition of the cross-correlation adopted in Eq. ([2]) differs from that found in other references (see e.g. 
Wapenaarl . 12004 . For deterministic or stationary random signals, the change of variable t 1 — > t' -\-t allows 



us to recover the form of cross-correlation used in other references. The difference is therefore purely 
notational. Since calculations are much easier to perform in the frequency domain, we will develop a 
theory for the cross-spectral density: 

/•OO 

C 13 {x b -x a -lo) = «u i (xfl;t>i(x il ;<' -t)*))e iut dt (3) 



In equation ([3]), the left-hand side does not depend on the variable t' for a stationary signal. Upon 
introducing the spectral representation of the random process Ui (x, t) : 



1 

Ui(x;t) = — / u^^e^du, (4) 

we obtain the following basic relation: 

(K(x B ;w')uj(x A ;w)*)) ^ 2nC tJ {x B ,x A ;uj)S{uj - uj'). (5) 

The delta function condition which appears in equation Q is characteristic of random processes that are 
stationary in time. We note that the usual Fourier transforms i(i(x; oj) may not be well defined since we are 
interested in stationary random fields. A physically satisfying solution is to define the Fourier transforms 
for a finite observation time T, and to take the limit T — > oo after averaging over noise sources. Instead, 
we make use of the probabilistic interpretation of the spectral representation (j4|), where the coefficients 
of the exponential is it self a random process, and the eq uality is to be understood in the probabilistic 
(mean-squared) sense ( Rvtov et al , 19891 Yaglom, 2004 ). Since only the cross-spectral density of the 



field is needed in this work, and this quantity can be defined by traditional Fourier analysis, the formal 
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representation Q suffices. A link between the physical and mathematical approaches is provided by the 
periodogram: 

Cij(x B ,XA;oj) = lim —({uj (x B ;u))uj(xA;u)*)}, (6) 

where uf(x; ui) represents the Fourier transform of the original signal observed over a finite time window 
T . As is well-known, Equation ([5]) expresses the fact that the periodogram is an asymptotically unbiased 
estimator of the cross-spectral density. 

The paper is organized as follows. In section [2] we introduce some important notations and present 
a simple derivation of Green's function reconstruction for a homogeneous medium illuminated by a set 
of randomly oriented and spatially uncorrelated forces. In section [3] we introduce the Dirac calculus 
which will be used to generalize our results to an inhomogeneous medium. The formalism is applied 
to demonstrate the reconstruction of Green's function of a general -homogeneous or inhomogeneous- 
elastic medium illuminated by a homogeneous distribution of randomly oriented forces, or by a wavefield 
at equipartition. In section 01 we introduce the T-operator for elastic waves and discuss its symmetry 
properties. In section[5j we explore the reconstruction of Green's function for an arbitrary inhomogeneity 
embedded in a homogeneous medium and establish several forms of the generalized optical theorem for 
elastic waves. In section [5J we establish the generalized optical theorem from the governing equations 
of scattering in the limit of infinitesimal absorption. In section [7j we consider a simple configuration 
composed of a single point scatterer and demonstrate the Green's function reconstruction including all 
the recurrent scattering loops. In section [51 our results are discussed in the framework of equipartition 
theory and compared to other works. 



2 Green's function reconstruction in a homogeneous medium 

In this section, we introduce important notations and present in a simplified context our approach to 
Green's function reconstruction. Although it is not the main theme of the paper, we give a short deriva- 
tion of Green's function reconstruction for elastic waves excited by rando m uncorrelated noise sources in 
a hom ogeneous medium. Our assumptions differ from the one adopted bv lSanchez-Sesma and Campillo 
(2006) who considered a set of uncorrelated plane P and S waves at equipartition. The complete equiva- 
lence between the two models of noise source will be demonstrated later in the paper. Let us first rewrite 
the correlation function ([3]) in terms of Green's function. The elastic wavefields recorded at two points 
and xb due to a distribution of forces can be expressed as: 



Ui(x B ;u}') = - JJJ G° fe (x B ,x ;w')A(xo;w')d 3 2o, (7a) 
Uj (x A ;ui)* =- III G° i (x j4 ,xi;w)*/i(xi;tj)*d 3 2; 1 , (7b) 



where G°(uj) is Green's function of a homogeneous, slightly dissipative elastic medium at angular fre- 
quency uj. Using the summation convention over repeated indices, G°(uj) is recognized as the fundamental 
solution to the equation of elastodynamics: 

(nodxidxiSij + (A + ^o)d Xi d Xj ) G° k (x 7 x';uj) + p (u> + i/2T) 2 G° ik (x 7 x';uj) = 5 ik 5{x,x') (8) 

and has the following analytical form: 

G° fc (x, x'; uj) = (g s (r; uj) + g n (r; uj)) (6 ik - hh) + (g p (r; u) - 2g n (r; «)) hh, (9) 

where Green's function has been split into transverse and longitudinal parts. In equation (JSJ we have 
introduced the notation r = x — x', the unit vector f in the direction of r, and the functions g p , g s and 
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^ik s r 



g a (T;u)=-- A j ( 10a ) 



gv(r;w)=-- - (10b) 



i / alp /D^kpV h, „ik s v i 



5 «(r; W ) =— 1- ^ + (e^ - e* k » r ) . (10c) 

The following short-hand notations have been used: u + = lu + i/^r, k p = uj + /c p and k s = uj + /c s , where 
c p and c s are the P and £ wave velocities in the background medium with Lame parameters Ao and /io- 
The definition of Green's function adopted in equation @ obeys an opposite sign convention to that 
usually found in the seismological literature. This choice will be clarified in section [3J In equation ©, 
we have explicitly separated the far-field terms g s and g p from the near field term g n . 

Equation (0 allows us to express the cross-correlation of wavefields measured at xa and xg in the 
frequency domain as follows: 

Gik{ x B,^o]^')G%{x A ,x 1 :ij)*((f k (x a ;w')fi(x 1 ;u)*))d 3 xod 3 x 1 (11) 

For spatially uncorrelated and randomly oriented forces, one further assumes that 

((/ k (x ; w)*}) = 2^4^(xo - Xl )S( W )(S( W - w'), (12) 

where S(u>) is the power spectral density of the forces. Therefore the key integral to be computed writes 

C , <i (xfl,x il ;«) = S(w) fff G° fc (x s ,x;w)G° (x A ,x;cu)*d 3 x. (13) 



— oo 



In coordinate space, Green's function has the complicated form @. However, using Fourier transforms, 
it is found that Green's function has the following exact expression in wavenumber space: 



G 



(k',k;a;) = g l (k;u)6(k' - k)fc. ( fc 3 + ff *(k;w)5(k' - k) - kkj) > ( 14 ) 



where we have introduced the following notations: 

ff'(k; w) = — 7— j 21, 2 i ■ 7~\ ' 5*(k;w)=— — 1 . (15) 

and the unit vector k is in the direction of k. In equation (|14[) . Green's function has been separated into 
longitudinal and transverse parts. Upon inverse Fourier transformation, both contribute to the near-field 
term in the spatial domain. In wavenumber space, the key integral is defined as 

3y(k / ,k;w) = ^4 III C iJ (xB,x^; W )e- jk ' XB+jkx ^ 3 x A d 3 XB. (16) 



Inserting the spectral representation of Green's function in equation (|13p . one obtains the following key 
integral in the wavenumber domain: 

p p p+oc 

C'y(k / ,k;w) = 5(w) /// G° ik (k',ko;u)G% (k, k ; W )*d 3 fc . (17) 



G 



Reporting expression (|14p into equation (IT71) . one obtains: 



Gyflc'.lqw) 



S(uj)t 
p oj 



lmg l (k)6(k' -k)kik 



S(uj)t 
p uj 



Imp*(k)5(k' ~k){5 lJ - kkj) 



(18) 



which establishes the reconstruction of Green's function of a homogeneous, slightly dissipative elastic 
medium. More precisely, equations ([3]), (fTTj) and (|18p allow us to obtain the following closed form 
expression of the cross-correlation function of two wavefields in the space-time domain: 

d 1 f T/2 

— lim — / ((ui(xB,t)Uj(xA,t -t)))dt' = 

at T—>oo 1 J 



2 Po (2tt) J_ 



duj 



e-^G?, (sbb.xa.w) ~ e^^C?, (x B ,x A ,u)} S (w) . (19) 



Hence, we retrieve the classical result that the temporal derivative of the cross-correlation of two wavefields 
is proportional to the difference of the retarded and advanced Green's functions filtered by the source 
power spectral density. Strictly speaking, the limit and integral signs on the left-hand side of equation 
(1191) arc superfluous. However, introducing the temporal averaging process makes a closer link with the 
measurement procedure. It is worth noting that the result (|19[) is valid under the assumption of equal 
absorption times for P and S waves. We will come back to this point in section [8] 



3 Green's function reconstruction in an inhomogeneous medium 
3.1 Formalism and notations 

To perform the calculations required for the proof of the elastic Green's function reconstruction in an 
inhomogeneous medium, it is convenient to introduce the Dirac formalism. This approach allows us 
to switch easily from coordinate to wavenumber space and is extremely efficient to perform compact 
formal derivations. In this paper, we work in a linear space of vector wavefunctions in which various 
representations can be introduced. It is convenient to think of this vector space as the tensor product 
of a 3-dimensional polarization space, and the usual space of scalar wavefunctions. In what follows, the 
symbol (u|v) will be employed to denote the scalar product between two vectors. For instance, the space 
of vector wavefunctions is equipped with the following scalar product: 

P P r + CO 

(u|v) = JJJ u^fvi^x, (20) 

which is inherited from the tensor product structure. In what follows, we will switch among various 
representations by introducing generalized orthonormal bases which verify the continuum normalization. 
The three most important bases for this work are listed below. 

(1) First, we define: \i, x) = |§j) ® |x), where the symbol ® denotes a tensor or direct product between 
vectors of the polarization and wavefunction space, respectively. The bra-ke t notations and definitions 
adopted here follow the treatment of vector fields in quantum mechanics (see Messiah . 19991 p. 549-551). 



The vectors e^) form an orthonormal basis of the 3-D polarization space. The simplest choice is to 
take |§i) as unit vectors along the axes of a right-handed orthogonal reference frame. We adopt the 
convention that unit vectors are denoted by a hat. The symbol |x) denotes generalized eigenvectors of 
the position operator. By forming the tensor product |e$) £3> |x), one obtains a wavefunction which is 
perfectly localized in space at position x and perfectly polarized along direction e^. The completeness of 
this set of orthonormal vectors is intuitively clear and the following normalization of the basis vectors is 
easily verified: (x, i\j, x') = (ej|ej)(x|x') = <5y <5(x — x'). Using this basis, we recover the usual definition 
of the recorded elastic wavefield in cartesian coordinates 

p p P+OO 

l u > = E/// Ui(x)\i,x)d 3 x. (21) 



7 




Figure 2: Conventions adopted for the definition of the elastic wave polarizations {Si, S2, P}- (ei, §2, 63) 
denotes a right-handed cartesian system, to which the local right-handed spherical system (6, <f>, f ) is 
referred to. An inversion of the propagation direction f — > — f, implies {6,4>) — > (#,—</>). The three 
polarizations Si, S2 and P are shown for the propagation direction f. 
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(2) Next, we introduce the symbol |i,k) = |e,) ® |k) which denote a basis of eigenvectors of the vector 
Laplace operator with eigenvalue — k 2 . The symbol |k) represents a plane wave with wavenumber k 
which obeys the normalization of the continuum. By forming the tensor product |e,) ® |k), one obtains 
a plane wave perfectly polarized along direction |e,-) and propagating in direction k. In the coordinate 
representation (1) introduced above, the eigenfunctions |z,k) write 

<M;,k> = |g^. (22) 

(3) Finally, we denote by |o;,k) the eigenvectors of the elastodynamic operator in free space L°, which is 
defined as 

(x,i|£°|u) = -(Aq + ^<d)d Xi d Xj Uj{x) - ^ dx J dx j u l {x). (23) 

The notational conventions adopted in Eq. (|23|) are as follows: the operator L° acts on the vector 
wavefunction |u) from the left to give a new wavefunction |v) = L°|u), and the right-hand side gives the 
usual expression of the wavefunction |v) in coordinate space. The minus sign adopted in definition (|23[) 
makes L° a positive operator, i.e., such that (u\L°\ u) > 0. This property is easily pr oved since (u|L°|w)/4 
represents the deformation energy of the solid fsee iBen-Menahem and Singh ( 1998 ) p. 31, for details). In 



coordinate space, the eigenvectors of L° are expressed as: 



Act ~ik-: 



where p Q denotes the polarization vector and k is the wavevector. Equation (|24p introduces the com- 
plete orthonormal set of linearly polarized plane P and S waves. In what follows, latin subscripts refer 
to cartesian coordinates (representations 1 and 2), while greek letters will be used to distinguish the 
contribution of the 3 different polarizations {Si, S2, P} in representation 3. The p" form a right-handed 
basis in polarization space which coincides with the usual spherical coordinate frame (6, c/>, f = k), where 
and 4> point in the direction of increasing 9 and 4>, respectively. Our conventions for the definition of 
the polarizations {Si, S2, P} are depicted in Figure [5J It is important to keep in mind that polarization 
vectors rotate with the incident wavevector k. Note that with our notational conventions, there should 
be no risk of confusion between representations (2) and (3). 



3.2 Resolvent formula and Green's function retrieval 

We now generalize the previous derivation to the case of an inhomogeneous elastic medium. We will 
consider simultaneously the case of a wavefield at equipartition, or a system driven by a white noise 
distribution of random forces. Our formalism allows us to treat the two cases on the same footing. Some 
basic facts about resolvents and their connection with Green's functions are recalled hereafter. Let us 
first introduce the operator L which pertains to elastic waves in inhomogeneous media: 

L = L° + V{ul), (25) 

where L° is the elastodynamic operator in free space introduced in equation ([23]) and V(u>q) is the 
scattering potential at the frequency ujq: 

MV{ul)\l^) = -5p(x)u$8(x-x.%i-d Xj SC ijkl (yL)d x J(x-x'), (26) 

where the usual summation convention over repeated index is used. Equation (|26l) defines the matrix 
elements of the scattering potential operator in the coordinate representation. This representation is 
most usually encountered in the Green's function retrieval literature. For simplicity, we assume that the 
scattering potential vanishes outside a bounded region of R 3 . To demonstrate that V(u>q) is Hermitian, 
we let V(luq) act on test functions f and g. Since the Hermitian symmetry of the term involving the 
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density perturbation bp is easily verified, we only consider the SCijki term: 
({\V(uj 2 )\g) = ///_ h^)*d X] (5C l3kl (*)dxMx))d 3 x 

(SCijki (x)d Xj /i(x)) * d Xk gi(x)d 3 x 



(27) 

(d Xk {SC ijk i (x)d Xj /i(x)) * gi{x)d 3 x 

(d Xk (5C kll3 (*)d Xj f 9l (x)d 3 x = (g\V(uj 2 )\f)*. 

In the derivation of equation (|27[) . we have used integration by parts twice and the usual symmetries of the 
elasticity tensor: Ciju — Ckiij = Cjiki = Cyifc. In the following, we assume that all Hermitian operators 
we have to deal with have unique self-adjoint extensions. Because the elastic tensor and density fields 
are real, relation (|27p implies that the matrix elements of the scattering potential in the representation 
(1) obey the following symmetry relation: 

<x,i|V(wg)|j,x / > = ^(x,x';w 2 ) = V^xjwg) (28) 

Let us introduce the resolvent of a linear operator L as follows: 

(XX — L)R(X) = X, (29) 

where X is the identity operator in the space of vector wavefunctions, and ImA ^ 0. Comparison of 
equations (|29|) and (O reveals that G°(ui) is intimately related to i?°(A), the resolvent of LP. By making 
the following substitutions A — >• po(uj + i/2r) 2 , L —> L , R — > G°(ui) in equation (|29|) . we can formally 
define Green's function at angular frequency uj with a small, finite absorption time r. This justifies a 
posteriori the definition of Green's function adopted in section [5J With this choice, G°(uj) is related to 
the resolvent of a positive operator with continuous spectrum on the positive real axis (see Equation (I23p 
and the discussion that follows). Let us now introduce in a similar fashion the resolvent and Green's 
function of the equation of elastodynamics in an inhomogeneous medium as follows: 

[XX-(L° + V(uj 2 ))]R(X)=X, (30) 
[po(w + */2r) 2 Z - (L° + V(lu 2 ))] G(u) - X. (31) 



An important point to be noted in the definition (1311) is the fact that w, the real part of the complex 
frequency, does not necessarily equal uq, the angular frequency at which the scattering potential is 
evaluated. For elastic waves, the only physically accessible Green's function is obtained when we set 
uj = ujq in equation (I31|) . Taking the matrix elements of equation (|31[) in representation (1), we recover 
the fact that Gj m (x, x'jwo) is the solution of Equation (|T|) with the force density fi — —Si m 5(x — x'), as 
it should. In equation (|3~Tj) . everything happens as if the scattering potential has been "frozen" at the 
frequency ujq. This mathematical trick allows us to use the standard tools of potential scattering theory 
an d facilitates interm e diate calculations. In elastodynamics, this idea has previously been introduced 
by iBudreck and Rose! (1991)- Eventually, all key formulas will be obtained for uj = ujq and will involve 



physically accessible quantities only. 

Let us consider the clastic wavefield at frequency uj\ excited by a distribution of forces and recorded 
at xs in an inhomogeneous medium: 



Uj(x B ; = - ^ J J J (x B ,i\G(uj 1 )\k,x 1 )(k,x 1 \f}d 3 x 1 , (32a) 
Ui(x B ;«i)=-5^ fff (p B ,i\G{ui)\P,Vi){flM\f)#ku (32b) 
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In equation (I32a[) . we have used the coordinate representation to express the field generated by the 
force distribution. But any other complete set of vectors can be inserted between G and /. As an 
example, in equation (|32b[) . we make use of this freedom by introducing a decomposition of the identity 
operator using the complete orthonormal set (3), which results in a hybrid representation. Indeed, 
the symbol (xb, i\G(tJi)\j3, ki) = Gip(xs, ki; cji) has its left foot in representation (1) and its right 
foot in representation (3). The quantity G^x^, ki; uj\) can be calculated from the usual coordinate 
representation of Green's function Gjj(x B ,Xi) by Fourier transformation over the variable Xi and by 
expanding the standard basis \ej) over the rotating polarization basis p^}. Clearly, the quantity /^(ki) = 
(/3,ki|/) gives the amplitude of the plane wave with polarization j3 and wavenumber ki, which is excited 
by the force distribution. As previously noted, if the force distribution is statistically homogeneous, the 
wavenumber decomposition of the field has a formal c haracter. Standard Fourier analysis can nevertheless 
be used after ensemble averaging ( Rvtov et al. ( 19891 )). Using equations (|32al) - (|32b|) . the cross-correlation 



between the waveficlds at x^ and xg can be written as: 

((uj(x B ;u;i)uj(xA;wo)*)) = 

[[[ (x B ^|G( Wl )|fc,x 1 )(xo,/|G( Wo ) t U,x A )(((/ fc (x 1 ; Wl )/ ; (xo; W o)*)))d 3 xod 3 a; 1 (33a) 



((ui (x B ; oji )uj (xa ; loq ) * )) 



E 

a,/3 



(x s , l |G( W i)|/3,k 1 )(ko,a!G(wo) t |j\x A )(((/0(k 1 ; Wl )/ a (k o ; W o)*)))d 3 fcod 3 fci, (33b) 



where G(wo)' denotes the Hermitian adjoint of G(wo). We now consider two possible physical assumptions 
regarding the excitations that drive the system. As in equation (TT21 , we may assume that the forces are 
randomly oriented and spatially delta-correlated. Alternatively, we may consider that the medium is 
illuminated by a set of plane P and S waves whose amplitudes constitute a white-noise in modal space: 

<(//j(k i; a;i)/ a (]co;wo)*» = 27n5 a/3 5(ko - ki)S(wo)*(wo - Wl ) (34) 

Assumption (|34l) may be interpreted as a requirement of equipartition: the system is illuminated by 
uncorrelated plane P and 5* waves coming from all possible directions with equal weights. This point will 
be further substantiated in section[S] Under which conditions such an equipartition state can be achieved 
in practice remains to be elucidated. Reporting Equations (TT2|) and (l34l) in Equations (|33al) and (I33bl) 
respectively, we find: 



Cij(x B ;x A ;u) ) =S(u} )^2 /// ( x b, i\G(u )\k, x)(fc, x|G(w ) t \j, xL A )d 3 x (35a) 

. J J J — OO 

Gjj(x B ;xA;wo) =5(a;o)E /// (x B ,i\G(ujo)\a,k )(ko,a\G(ujo)''\j 7 x A }d 3 k (35b) 



=S , (w )(x s ,i|G(o;o)G(a;o) t |i,x J 4), (35c) 

where in the last equation, we have used the completeness of the bases (1) and (3), respectively. In both 
cases -white noise distribution of forces, or wavefield at equipartition- we have reduced the computation 
of the cross-correlation of two random wavefields to the evaluation of the matrix elements of the operator 
G(cjo)G(cjo)^. In fact, from the defining property of Green's function as the resolvent of a self-adjoint 
operator -equation (13"TT) -. the following identity follows: 

GMGM 1 =R(po(u + i/2t) 2 )R( Pq (lu - i/2t) 2 ) 

= - — (R{po{ojo + i/2r) 2 ) - R{p {oj - «/2t) 2 )) ( o 6) 

= --^(GM-GWt), 
Zipoujo 
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with R(X) the resolvent operator introduced in Equation (|30p . The key point in the above derivation 
is the use of a fundamental operator equation known as the "first resolvent formula" in the second line 
of Equation (1361) . The resolvent formula is valid for a broad class of ope r ators (closed linear operators) 
including the differential operators considered in this work ( Richtmverl Il978l ). The last step follows 
simply from the Hermitian character of the elastodynamic operator L° + V(u>q). Taking into account 
the power spectral density of the source and calculating the matrix elements of the right-hand side of 
equation (f36]) yields: 



d 1 f T/2 

-r. lim - / ((u i (x B ,t)u j (x A ,t -t)))dt' 
at l — >oo l j _ T 12 



2 Po {2n) J_ 



du 



e-^Gij (x b ,x a ,oj) - e-^-QGjiix^XB,!*)] S (w) , (37) 



where the Hermitian symmetry of G(uS) has been used. When the reciprocity theorem applies, the 
integrand in the right-hand side of equation (|37l) is recognized as the usual imaginary part of Green's 
function. Thus, the derivative of the cross-correlation tensor of random signals yields a filtered version 
of the difference between advanced and retarded Green's functions. Since G can be substituted with the 
Green's function of any self-adjoint operator, this implies the reconstruction of Green's function for other 
systems governed by linear wave equations. Derivations of Green's function reconstruction based on a 



differe nt operator formalism have been further developped by IWapenaar et ali (|2006l ) and ISnieder et al 



(2007). Our derivation highlights the fact that a wavefield at equipartition and a random distribution 
of white noise sources are indistinguishable from the point of view of cross-correlation functions. The 
complete equivalence between the two types of excitations will be demonstrated in section |8l The general 
result (f36|) does not constitute, however, the end of our investigations. We will show that interesting 
information on the scattering of elastic waves can be obtained by expressing the Green's function of an 
inhomogeneous medium using the formalism of scattering theory. To do so, basic notions of scattering of 
elastic waves are recalled in the next section. 



4 Summary of scattering of elastic waves 



Our starting point is the observation that Equation (|31[) is formally identical to a Schrodinger equation, 
with a scattering potential V(u>q). This allows us to use all the arsenal of scattering theory. In particular 
we can introduce the retarded Lippman-Schwinger eigenvectors of th e operator (L° + V (un)) which satisfy 



the following equation in the infinitesimal absorption limit t — > oo ( Economou ( 20061 )) 



|?/; a (fck)} = \a, fck) + G°(c a k)V(L>Z)\iP a (kk)), (38) 



with fc > 0. It is worth noting that in equation (|38j). the potential is "frozen" at the frequency luq, and 
is scanned by plane P and S waves with all possible wavevectors fck with fc > 0, and eigenfrequency 
u = c a k > 0. These eigenvectors do not depend on the sign of loq. For elastic waves, the physical 
wavefunction is obtained when the condition lj = ujq — c a k is satisfied. For later developments it is 
nevertheless convenient to introduce the more general Lippmann-Schwinger eigenvectors p8[) . All the 
scattering properties of an inhomogeneous elastic medium are conveniently encapsulated in an operator 
denoted by T^o) which is defined as follows: 

T("°)(u;) = V(^ ) + V{ul)G{u)V{ul). (39) 

In equation (f3"9"|) . we have carefully distinguished between the angular frequency of the target -uj q - and 
of Green's function -u>-. Clearly, the physical result is obtained when uj equals luq and we introduce a 
special notation for this case: 

T° = T^o\lo ) (40) 
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For simplicity, we have dropped the frequency dependence in this definition. The superscript reminds us 
of the fact that 10 = u)o- The definition of the T^o)_operator adopted in equation (|39p is in keeping with 
our definition of Green's function of an inhomogeneous medium (see equation (|3ip ) and will facilitate 
intermediate calculations in section [5] Using the T^o) -operator, we can rewrite the Lippman-Schwinger 
equation as follows (r — s> oo, k > 0): 

\ip a (kk)) = \a,kk) + G°(c a k)T^« ) (c a k)\a,kk), (41) 

which further implies: 

V(u^)\ip a (Jtk)) =T^?>\c a k)\a,kk) (42) 

Green's function of the heterogeneous medium can also conveniently be rewritten in terms of the T^°°^- 
operator and the free space Green's function alone: 

G(u) = G°(w) + G°(w)tK 2 )( w )g°(w). (43) 

Since our approach is based on an analogy with quantum scattering, a word of comment about the 
sign of wo is in order at this point. The problem of negative frequencies is easily solved because we 
consider real fields in the space-time domain. This imposes the following Hermitian symmetry on the 
matrix elements of G and T in coordinate representation: 

G^(x,x'; -wo) =Gy(x,x'; w )* T^W'; -w ) =T^ (x, x'; w )* (44) 

As a consequence, all results obtained for positive frequency are readily translated to negative frequency 
using the symmetry conditions (|44[). 

In the far-field of the scatterer, i.e. in the limit x 3> D with x the observation point and D the 
dimension of the scatterer (Figure [1]), we may approximate the free space Green's function as follows: 

G° (x, x') ss — (Sij — x t Xj) - x^j, (45) 

J inpocfx 47rp CpX 

where x' denotes an arbitrary point of the scatterer. The position vectors are measured from an arbitrary 
origin located inside the scattering region and x is a unit vector in the direction of x. The asymptotic 
formula (|45j) allows us to express the Lippman-Schwinger wavefunction in the usual coordinate represen- 
tation as follows: 

<i,x|^k fe )) = (Ve^- + grffelM ffe*"-) , (46) 

valid in the limit x — > oo, where x = |x|. In equation (|46j) . we have introduced the scattering amplitudes: 
r^(k Q , k b ) = - (2 ^\ 2 (a, fc Q k a |T^ 2 >M|/3, k^ b ) 

o a- r (47) 



47rpo(c a ) 5 



and the notations: k a ^ — u>o/c a _p > 0. In the last line of equation (14"T1) . we have introduced "on-shell" 
matrix elements of the T°-operator. In general, the matrix elements (a, k Q |T^o' (w)|/3 ,k;,) are said to be 
"on-shell" when the condition uS 2 = (? a k 2 a — c^k 2 is satisfied. Although this property does not depend on 
the fact that we are on the shell ui = ujq, we eventually specialize our results to this particular case. The 
scattering amplitudes contain information on the scattering pattern for all possible mode conversions. For 
instance, taking /3 = 1, a = 3, corresponds to the physical situation where a linearly polarized incident 
shear wave is scattered into a compressional wave. Note that the definition is valid for an arbitrary 
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anisotropic inhomogeneity. From the symmetries of the scattering potential: (x, x'; Wq) = V,i(x', x; Wg) 
and the reciprocity of the free space Green's function G^(x, x';u;) = G^(x', x; cl>), we deduce a similar 

reciprocity relation for the T^ w ^ operator: 

^ 2) (x,x'; W )=r 7 ( r° ) (x',x; W ). (48) 

Using equation (1481) , our formalism allows us to recover well-known reciprocity relations for the scattering 
amplitudes: 

T { a f(K,k b ;u) =E(i / // J?e- ik ^T^\^y-,Lj)q^ k ^'d 3 xd 3 x' 



1,3 



y- 



(2tt) 3 



q?e ik >-*'T!>f°\x\x;u)p?e- ik °-*d 3 xd 3 x' ( 49 ) 



3 / / /_„ ^3 3 l 



=sgn(a,/3)T£ o) (-k b! -k Q ; W ), 

where sgn(a,/3) = —1 for (a,/3) 6 {(1, 2), (2, 1), (1, 3), (3, 1)} and +1 otherwise. In equation (j4*9"|). the 
matrix elements are generally off-shell because the condition uj 2 = c„fc^ = needs not be satisfied. 
This relation is in fact far more general than what is required for our purposes. We will therefore 
specialize it to the on-shell case with the condition u = wo- Note that the sign change -sgn(a, j3i)- is a 
mere consequence of the convention adopted to define the polarization vectors and illustrated in Figure 
[2j Indeed, upon a change of the wavevector k a — > — k a , the polarization vectors transform as follows: 
(p^P^P 3 ) (p 1 ; — P 2 , — P 3 )- As an example, let us consider an incident S wave with polarization 
vector q 1 and a scattered P wave with polarization vector p 3 . Application of the symmetry relation (|4"9"f 
to equation (|47[) yields the following reciprocity relation between the scattering amplitudes: 

C 2 / P - Sl fc 6 ) = -C 2 J S ^ P (-k b , -ka), (50) 



a relation which has been noted by a number of authors in the past ( Varatharaiulu , 1977 ; Dassios et ali 



1987tlAkil . 119921 ). Interestingly, this reciprocity relation shows up without any effort from the formalism we 
have introduced. Some physical implications of relation (|50[) will be discussed in section |8l in connection 
with the role of equipartition in Green's function retrieval. With the help of our scattering formalism, 
we will explore the consequences of Green's function reconstruction in a scattering medium in the next 
section. 



5 Green's function reconstruction in a scattering medium 

As noted in section ([2]), Green's function reconstruction from homogeneous random sources is a con- 
sequence of a basic operator identity. We now explore the consequences of this result in the case of a 
scattering medium. The identity satisfied by Green's function of the heterogeneous medium writes: 

GMGM 1 = (G(wo) - CM 1 ). (51) 

Using equations (|43|) we can express the left-hand side (L.H.S.) of equation (f5Tj) as follows: 

GMGK)t =GV )GV )t + G ( Wo )G ( Wo )tr°tG ( Wo )t (52) 

+ G MT°G MG°M t (53) 
+ G°(io a )T G a (uj )G a (Lu o yT^G a (io a y. (54) 

Similarly, we obtain for the right-hand side (R.H.S) of equation (|5T|) : 

G(«o) - G(w )t = G°M - G°(o;o) t + G°(u )T G (lu ) - G°(wo) t T ot G (w ) t . (55) 
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Upon inserting relation (j3"o| into equation (|5"i|) . Green's function reconstruction formula ([5T|) yields the 
following identity for the T operator: 

I (T° - T°t) = -^T°G (lj )G°(uj ^T^ (56) 

In order to elucidate the physical meaning of equation (|56[) . let us take matrix elements in the represen- 
tation (3). We find: 

l(a,k | (T°-T°t) |/?,k b ) = | (T^(k a ,k fc ) -T° Q (k 6 ,k a )*) , (57) 

on the L.H.S., and: 



(a,k a |T°G°(a;o)G ( W o) t T ot |/3,k 6 ) 



wo. /•/•/■ (a,k a |T | 7 ,k)( 7 ,k|r°t|^k b ) d 3 /; 



on the R.H.S., where we have made use of the spectral decomposition of G°(u>o), valid irrespective of the 
sign of luq : 

= £ /// , J 7 '^'- (59) 
^ J J JR3 M^o ~ S fc + iW o/ r ) 

On the right-hand side of equation (|59p. the symbol |7,k)(7,k| is an operator which projects the wave- 
functions on the subspace of plane P and S waves with wavevector k. We will now turn equation (|58p into 
a more physically appealing form. First, we note that the representation of the product of Green's func- 
tion G° (u>q)G° (wo)^ is left unchanged by the substitution k — > — k. Next, we make use of the reciprocity 
relation (|49|) to rearrange the numerator as follows: 

(a, k a |T°| 7 , -k)< 7 , -k|T°t|/3, k 6 ) = sgn(a, (3)(f3, -k 6 |T°t| 7 , k)( 7 , k\T°\a, -k a ), (60) 

where we have used sgn(a, 7)sgn(7, f3) — sgn(a,/3). The following substitutions: k a — >• — k b , k b — > — k a , 
a —>/?,/? —> a, leaves equation (|57|) unchanged except for an additional factor sgn(a, (3) which cancels 
the same factor on the right-hand side of equation (|60p . We therefore arrive at the following general 
identity: 

which we will call the generalized optical theorem in the presence of small finite absorption. Equation 
(foTjl is a general identity which must be obeyed by the "off-shell" matrix elements of the T°-operator 
for a general scatterer in an elastic medium. The matrix elements are said to be off-shell because k, 
k a and k(, do not have to verify the shell condition ujq = c 2 k 2 . The physical meaning of these off-shell 
contributions will be discussed below. Equation (I5T|) is a central result of this paper, from which all forms 
of optical theorem derived in the literature can be deduced. The justification of the term "generalized 
optical theorem" is made clear when we consider the limit of infinitesimal absorption r — > oo. Using 
standard results from distribution theory, we find: 



Wo /(por) 



cU 2 ) 2 + uj 2 /t 2 2p w c 7 . 



(62) 



Upon inserting the previous relation in equation (|61|) and integrating over the wave number fc, we arrive 
at: 

\i C**»"»> -nJW) = -^^Jl I* *) C (^t, - rt, (63) 
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which is the form of the generalized optical theorem in a medium with infinitesimal absorption. In 
Equation (|63l) . the symbol cPk denotes the element of solid angle in wavenumber space. Again, off-shell 
terms of the T matrix appear on both side of equation (|63[) . On the right hand side, two arguments are 
on the shell c 2 fc 2 = uIq. Physically, off-shell terms are related to the physical situation where source and 
detection take place in the near field of the scatterer riPeierlslll979h . Thus the reconstruction of Green's 



function for arbitrary locations of the receivers in the medium imposes symmetry relations that pertain 
to the off-shell T-matrix. More examples of the optical theorem for electromagnetic and quantum waves 
can be found in the book of lNewtonl <|2002h . If we now impose a = f3 and further restrict all wavenumbers 
to the shell uiq, we arrive at classical forms of the optical theorem: 



ImT" 



2po 



r3 



t: 



7 a 



w . 



ui 0f 
— k a 

Cry 



d 2 k 



(64) 



As discussed in iPeierlsl (|1979l ) , the on-shell matrix elements describe the angular dependence of the 
scattered waves in the far-field of the heterogeneity. The physical meaning of equation for elastic 
waves will be further elucidated in section [6] It is important to keep in mind that all the relations 
obtained so far have been based on the Green's function reconstruction theorem. In the next section, we 
make the link with the scattering theory developed in section |4] 



6 Generalized optical theorem: formal derivations 



In this section, we provide derivations of the various forms of optical theorems obtained within the 
framework of Green's function reconstruction. First, we note that Equation (|64[) which is the on-shell 
version of equation (|63|) . can be obtained straightforwardly from the scattering theory developed in section 



1 

2i 



1 a8 



1 

2i 
1 

2i 



^k —k h 

C a CS 



1 

w . 



■T, 



1 

2i 



8a \ 

\CB C, 

(a,^k a |T°|/3,^k,)-(a,^k a |T°t|^^ k6 

C0 c a Cp 



Co \ Cf) 



c /3 



ir ( -k a ) \v{uD\^ ( ^k b )) - (r [ ^k«k a ) \v(ui)G°{^v{ui)\p, ^k b ) 

C a ) \C0 J \C a J C/3 



UJQ . 



1 

" 2i 



^k a |FK 2 )G ( w )nw 2 )|^Pk h 



(r (^k a j w^d (gV) - gvy) v{u>d\^ ^k b j) 



i 

2i 



(a, ^k a |T°t (G°( w °) - G°(c°)t) T°\f3, ^k b > 



TTUJq \ - 1 
9n„ 2-^1 77~ 



2p - ^ 7 



4tt 



C 



^(^^Wo.lSSfi,^) Sk (65) 



where the infinitesimal absorption limit (r — > oo) is understood throughout. The first equality follows 
from the definition of the T°-operator. In the second equality, we make use of the relation (l4"21) be- 
tween Lippman-Schwinger eigenvectors and plane waves. The third equality follows from the Lippman- 
Schwinger equation (|4Tj) . The fifth equality makes again use of (|42|) . We recover the generalized optical 
theorem (|63|) after inserting the spectral decomposition of the free space Green's function ([59)) in the 
sixth equality. Some on-shell forms of the optical theorem similar to (|64|) -(|65 j) have been previously 
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derived in Lu et al. ( 2011 ) for a homogeneous and isotropic heterogeneity, and Budreck and Rosel ( 1992f) 
for a general heterogeneity. These formulas are less general than the fully off-shell expressions (|6Tj) and 
(|63p . The physical meaning of the formula (|64l) can be further elucidated by introducing the scattering 
amplitudes defined in (|4"T)) . From the definition of the scattering cross-section as the ratio between the 
total energy scattered per unit time, normalized by the energy flux density of the incoming plane wave: 



(66) 



we arrive at the following simple formula: 

c Q (k a ) 



4tt 

ka 



imr^ Q (k a ,k a ), 



(67) 



which has been derived by other authors in the past ([Varatharaiulul ([19771 ); iDassios et ali (|1987j)). The 
physical meaning of equation (1641) is now clear: it expresses the conservation of energy in the scattering 
process. Energy conservation is realized thanks to the interference between the incident wave and the 
wave scattered exactly in the forward direction in the far-field. 

To prove more general versions of the optical theorem, we invoke the spectral decomposition of 
Green's function of the heterogeneous medium in terms of the Lippman-Schwinger eigenvectors. We 
make no attempt to justify rigorously this decomposition and proceed comp l etely formally. Mathematical 
founda t ions o f this approach can be found in the book of iReed and Simo 3 (|l979h for quantum waves, of 



Ramml (|1986l ) for classical waves. A generalized eigenfunction expansion in te rms of Lippman-Schwinge r 



eigenvectors of the elastodynamic equation may also be found in the work of iBudreck and Rose! (Il99lh . 
Making use of the eigenvectors of the operator L° + V(u!q) introduced in equation (|38|) . we expand Green's 
function G as follows: 



\^{kk)){ip" 1 {kk)\ 



k l dk. 



In the case of quantum waves, this decomposition is correct if the scattering potential does not create 
bound states. Using the formal definition of the T^ u o) -operator ([59")) and the spectral representation ([6"8")) . 
we find in the limit (r oo): 



i (T^(k a: k b )-T° a (k h ,k a y) = i (<a,fc a |U(w 2 )|/?,fc 6 ) - <a,fc a |V(w 2 )|/3,fc 6 }) + 

i ({a,k a \V{ul)G{uiv)V{ujl)\p,k b ) - (a,fe a |y( W 2 )G(wo) t F(w 2 )|/3,fc b }) 



POT ^ 



4tt 



d 2 k / fc 2 (a,fc a |U(w 2 ) 



\^{kk)){^{kk)\ 

[ (w 2_ c 2 fc 2 )2+w 2 /T2) ] 



V(uj 2 )\l3,k b )dk 



+ 00 



— y 



d 2 k 



c 7 w 



(5(k - w /c 7 ) + S(k + wo/c 7 )) 



x (a,fc a |T w o(c 7 fc)|7,fcfe)(7,fcfc|T""(c 7 fc)|^,fc b )dfc 



7TW X - 1 



c 7 

7 7 



4tt 



W 0f 



drklg} I h a , — k; c 7 fc = w ) Tp° [ k b , — k; c<k 



W 0f 



w 



TTUJo \ -> 1 



c 7 

7 7 



d 2 fcr 



w 



k, fe a 



J 7/3 



W 



k, feft 



(69) 
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which is identical to equation (1551) . In the derivation of equation (jSHJ), we made use of the limit (p32"j) 
and of the reciprocity relations for the T^o) -operator. In the third equality, we have made use of the 
Lippman-Schwinger eigenvectors introduced in equation (|41[) . It is to be noted that the previous results 
are correct irrespective of the sign of ujq, which is important when dealing with classical waves. 



7 Application to a point-scatterer model 

We now illustrate the process of Green's function reconstruction with a simple example of a 3-D het- 
erogeneous medium composed of a single delta-like scattering potential. To render multiple-scattering 
calculations easier, it is often convenie nt to introduce a point-scatterer model in the same spirit as in the 
scalar case dMargerin and Satol l201lh . Mathematically, point-scattering poses some serious difficulties 



( de Vries et all , 19981) . Physically, the point scatterer model is obtained by regularizing the Born series 



for a scattering potential of the form: 

V(w 2 H 7 (^|xo}<xo|, (70) 

4) 



where Xq is the location of the scattering center, and the physical significance of the parameter 7(0^ 



will be further discussed below. The scattering potential V(oJq) has to be understood as a tensor product 
between two operators. The symbol I denotes the identity operator in polarization space. It simply leaves 
unchanged the vectorial part of the wavefunction. The symbol |xo)(xo| is a projection operator which 
picks out the value of the wavefunction exactly at position Xo, hence the name "point-scatterer". We will 
find that the scattering potential V(u>q) represents adequately a small scatterer with perfectly correlated 
perturbations of the density and P and S wave velocities according to: Sc p /c p — Sc s /c s = —Sp/2p, 
(SX/X = Sfi/ fi = 0). No assumption on the smallness of the perturbation is made. In the representation 
(2), the scattering potential writes 



^•(k,k';a; 2 ) = ^ e < -( k '- k )%. (71) 



To evaluate the scattering properties of such a perturbation, we make use of the formal Born series for 
the T-operator of elastic waves: 

T° = V{uj 2 ) + V{u$)G Q (uo)V(u$) + V(Lu 2 )Go{uJo)V{u]l)G a (^)V{u]l) + ... (72) 

and evaluate the matrix elements in the representation (2) 



(k,i|T°|j, k') = ^(k,k'; Wo 2 )+^ y y y^(k 1 *|V(wg)|* lj k x ><t 1 ,ki|G^(w )V(a;g)|j,k / >d 3 Ai+ 

III <M|t>o)l*i,ki><ki,*il^(^)v^^ 



+ ■•• (73) 

The need for regularization stems from the fact that repeated visits of the same scatterer involve evalua- 
tions of the free space Green's function for coincident source and receiver, which clearly diverges. We now 
show that for a scattering potential of the form ([701) . the Born series (|T2"]) can be summed analytically 
after the return Green's function is made finite. To make sense of expression (|73p . we need to evaluate 
the matrix elements of the operator Gq(ujq)V(ujq): 



(74) 



<M|GVo)t'(wo)|j,k / ) =E /// <k,i|G°(wo)|ii,ki>(ii,ki|^(«g)|j,k'>d 8 fci 

J J J 

7(w 2 )e 4X °-( k '- k ) 




p (w 2 - (c p k) 2 + iujo/t) p (uj% - (c s k) 2 + iu /r) J ' 
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where we have made use of equations (TH| and (|TTj) . Reporting expression ([71)) in the first integral of 
Equation (|73[) we obtain: 



h = 



7(^o)e 



2\ ix -(k'-k) 



kiki 



T-kxki 



(2tt) 3 

7(wg)e JX °( k '- k ) X47r 7 (o;2) 



po(^o ~ (cp^i) 2 + iloo/t) po(wq - c 2 fc 2 + iwo/r) 



k\dki 



(2tt) 3 

where we have used: 



3(2tt) 



Po(wq - c 2 fc 2 + iuiq/t) 



k\dk\ 



Pq{ujI - c 2 s k 2 + iuJo/i 



kikid 2 fci = 

4tt 



Similarly, for the second integral of equation (|72|) . one obtains: 



/ 2 = 



7(^o) 
(2tt)3 ( 



ixo(k'-k) 



X47r 7 (q;g) 
3(2tt)3 



kfdki 



Poi^o ~c 2 v k\ +iu /t) 



kfdki 



A)( w o ~c 2 k 2 +iuj /t) 



(75) 



(76) 



(77) 



From equations (f75j) -(f77 |) . we infer that the Born series is in fact a simpl e geometric series. The last inte- 
grals i n ([75|) - d771) diverge and need to be regularized. We proceed as in Ivan Rossum and Nieuwenhuizen 
(1999), by introducing a large momentum cut-off to avoid divergence of the real part of Green's function. 
Let us illustrate this procedure for the P modes. We rewrite the first integral on the right-hand side of 
equation (|77|) as follows: 



J 



Poet 



k\dk\ 



k 2 k 2 -\- ikp/lp 



—J* 

PoCp J 



dk\ 



(kp k\ + ikp/lp) 



(78) 



kj 



Poc 2 , 



dki, 



where k p = ujq/c p and l p = c p t. Note that terms of order 1/Z 2 have been neglected in the numerator, 
which implies that absorption is sufficiently small. In the first integral, we have subtracted out the 
divergent part. In the second integral, we introduce a large momentum cut-off with a scale length a to 
obtain: 



(k p + i/2l p ) 2 dk 1 



l/a 



Pocl J k 2 ~k 2 + ikp/lp p cl 


poc 2 , \a Alp 2 



dki 



(79) 



where k p = ujq/c v . Combining equations (|75|) -(I77 |) and (|79|) . one may rewrite the formal series expression 
of the T°-matrix given by (|72j) in the form: 



T£(k,k') 



e ,x »( k ~ k )i(w )<5 ! . 



where the scalar i(wo)-rnatrix can be expressed as: 



t(uj ) 



(27T) 



1 -7(^)^0)' 



(80) 



(81) 
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and we have introduced an effective potential: 



Up ^0 
2 C 3 + cf 



(82) 



a \ ct ct 



where the length a stems from the large momentum cut-off. The result (|8Tj) is valid irrespective of the 
sign of wo- Equations (|80p -(f8T j) offer a generalization of the poin t-scattering m o del to vector elastic 
waves. The regularization procedure is not unique as illustrated by lde Vries et a l. (1998) for scalar and 



electromagnetic waves. The most important point is the preservation of the local imaginary part of the 
free space Green's function: 

+ (84) 

which is proportional to the density of states, i.e. the number of modes per unit volume and frequency. 
We note that the imaginary part of the free space Green's function at (x, x) is an equipartition mixture 
of P and S waves and does not depend on the absorption time r. The result (I81j) for the i-matrix is valid 
for a small, finite amount of uniform absorption inside the medium. Note that because the absorbing 
properties of the scatterer match those of the embedding matrix the scattering potential is real. 

To analyze the properties of our point scatterer, we make use of the relation (|47p providing the relation 
between the scattering amplitude and the T°-operator. The scattering pattern is obtained by sandwiching 
the T°-operator between two plane waves, i.e. we take the matrix elements of T° in representation (3). 

r^(k a ;k fc ) oc (k a ,a\T°\k b ,p) cx <p Q |I|cf ) = p a ■ , (85) 

where p Q ,q^ denote the polarization vector of the scattered and incident wave, respectively. Let us 
specialize the last result by introducing a spherical coordinate system where the polar angle is measured 
from the axis pointing in the direction of the incident wave with wavevector k^. For the different scattering 
type, we obtain the following angular dependence: (1) P to P: cos#; (2) P to SV and SV to P: sin#; 
(3) SV to SV: cos6> , (4) SH to SH: no angular dependence (in and out polarization vectors are 
parallel), where SV and SH denote transverse polarizations in and perpendicular to the scattering plane, 
respectively. All other coupling terms are zero. These results agree with the scattering pattern for a pure 
density perturbation (constant Lame parameters) in the low-frequency or Rayleigh regime. The scattering 
pattern is anisotropic but non-preferential, i.e., there is equal amount of forward and backward scattering. 
This property is typical of a scatterer whose size is much smaller than the wavelength. There is no mode 
conversion in the forward direction as required by symmetry consideration. Reporting expression (|80[) in 
the usual form of the optical theorem (|67p , we find the following relations for the scattering cross-sections 
in terms of the t-matrix: 

1 -Im— (86a) 



p c p ujq 

— Im— , (86b) 



which implies: 

Cp(T p = c s cr s . (87) 

The relation (187[) is valid for a point scatterer but not generally true. Finally, let us identify the parameter 
7(wq) with — SpujqV, with V the volume of the scatterer and 5p the density contrast between the scatterer 
and the matrix. By substituting expression (f8Tj) into expression (|86a|) . and taking the limit ojq — > 0, we 
find: 

- 2V fj_ 1_\ (5p^ 2 - 4 ~ 3 
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with r the scatterer radius. These results concur wit h previous analytical solutio ns of the scattering 
problem for a spherical inclusion in the Rayleigh regime ( Korneev and Johnson , 1993). We have therefore 



defined a physically meaningful point-scatterer for elastic waves. It is important to note that relation 
(|87p pertains to a localized perturbation of the density at fixed Lame parameters A and /z. 

The final step of our investigation of the point scatterer concerns the verification of the complete 
identity (|6Tj) . On the right-hand side, we find: 

2l {-Lapika^b) - lp a {k b ,k a ) ) = ^— ^ (89) 

The left-hand side can be calculated as follows: 
w ° V- fff ^(k,k b )T 7 ° Q (k,k a )* 



Por^JJJ (^-(c 7 fc)2) 2 + c 2 /r2 



Wo ^ fff |t(a;o)| 2 e < *°-( k »- k °> (p° | 7 , k) (7, kjq^) i3fc 



- Wo |^ )| 2 p"-q^ /l 



(27T) 3 67rp W 2 4 



(90) 



a result valid in both absorbing (finite t) and non-absorbing media (r — >■ 00) . Fr om equation ([84)1. we 
deduce a relation analog to the scalar case (jvan Rossum and Nieuwenhuizen , 1999() : 



Imi(a;o)J= |i(w )| 2 ImG°(x, x; w Q ). (91) 

A straightforward calculation shows that relation (1911) is indeed verified by the t-matrix obtained from the 
regularization procedure. With this simple example we have therefore illustrated the relation between the 
reconstruction of the elastic Green's function and the conservation of energy in scattering. These results 
can be generalized to more complex multiple-scattering media. Such developments will be presented 
elsewhere. 



8 Discussion and Conclusion 

The reconstruction of Green's function of a homogeneous e lastic medium from random noise sources was 
previously studied bv lSanchez-Sesma and Campillol (l2006h . These authors used a decomposition of the 



random wavefield into a sum of uncorrelated plane P and S waves at equipartition. In this paper, we 
gave a precise mathematical formulation of the equipartition condition in equation (1341) . Using the Dirac 
calculus, we can show the complete equivalence between the equipartition state and an excitation of the 
system by randomly oriented and homogeneously distributed uncorrelated forces. Noting the relations: 

/j(xi;wi)/ fc (xo;a;o)* =(Z,x 1 |/(wi)}(/(wo)|fe 1 xo) (92a) 
//KkijwO/aOcojwo)* =(^,ki|/(«i))</(wo)|a,ko>, (92b) 
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Figure 3: Energy balance of a random wavcfield in a homogeneous isotropic elastic medium containing a 
scatterer. The scatterer is illuminated by uncorrelated plane P and S waves coming from all possible space 
directions. The total energy carried by P waves increases because an S wave propagating in direction 
n' can be converted into a P wave propagating in direction n. The reciprocal situation shows a P wave 
propagating in direction — n and converted into an S wave propagating in direction — n'. This scattering 
process entails a decrease of the total P wave energy. For an incident wavefield at equipartition the two 
contributions cancel exactly. The same argument applies to the total S energy. 
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we obtain: 



«//j(ki;wi)/ (ko;a;o)*» = 

((((/3,ki|Z,xi)(?,x 1 |/(wi))(/(w )|fc,xo)(fc,x |a,k })))d 3 a;od 3 xi 



E 

k,l 



2ttS(ujq)S(u!o — u>i) 

^ / / / S(-k - xi)^ fc i((/3,ki|fc,x )(?,Xi|a,ko))d 3 a;od 3 xi 



x 



(93) 



=27t5 , (w )'5(wo - w i)E /// (/3: k i| fc , x o)(fc,x |a,k )d 3 a;o 

J J J 

=27r5'(wo)'5(wo - wi)<5 a /3(5(ko - ki) 

Equation (|93[) . demonstrates that from the point of view of correlations, a white noise distribution of 
random forces and a waveheld at equipartition are equivalent. We may also illustrate this point with 
physical arguments. First we note that the radiation of energy by a single force can be evaluated from the 
far-field solutions. For the P-wave radiation, we multiply r 2 pqu>qC p by the squared modulus of Green's 
function g p (r;uio) as given by Equation (jlObl) . The procedure is the same for the 5* wave radiation, 
except for the substitution c p — > c s . It follows that the ratio between the total S- and P- energy 
radiated by a force is 2(c p /c s ) 3 , which rigorously coincides with the ratio of the equipartition state. A 
homogeneous distribution of randomly oriented forces will therefore generate P and S waves propagating 
in all possible directions with a local energy density of P and S waves satisfying the equipartition 
relation E s /E p = 2c 3 /c 3 . Thus we find that the force distribution (fl~2"|) guarantees the existence of the 
equipartition state and the reconstruction of Green's function. It is important to note that the model 
of absorption that we have developed implies equality of the intrinsic quality factor between P and S 
waves. This assumption may not always be fulfilled in practice. If the absorption times of P and S waves 
are different, the P-to-S energy ratio differs from that of the equipartition state and Green's function 
reconstruction will be imperfect. 

The requirement of equipartition fo r the successfu l reco nstruction of Green's function of a heteroge- 



neous medium has been widely noted. iDerode et all (|2003l ) have studied numerically and theoretically 



the reconstruction of Green's function of a heterogeneous, open medium illuminated by either random or 
deterministic sources. They gave a criterion based on time-reversal arguments to discuss the optimal ar- 
rangement of the sources. The reconstruction of Green's function of a heterogeneous medium illuminated 
by an incident random w avefield at equipartition has b een considered theoretically bv I Weaver and Lobkis 



(2004) for scalar waves. Sanchez- Sesma et have demonstrated that when a cylindrical scat 

terer is illuminated by random plane P and S waves at equipartition in 2-D, the cross-correlation of the 
wavefield allows the reconstruction of the complete Green's function, including the waves scattered by 
the cylinder. In the case of a distribution of sources on a re mote surface enclosing the receivers and a 



3-D isotropic and homogeneous heterogeneity, iLu et all (|2011l ) showed that the wavefield incident on the 



scatterer must obey an equipartition relation. Their derivation supposes that the medium is free from 
absorption. In this work, we have generalized these results to an arbitrary, localized heterogeneity and 
have given an application to a point-scatterer. To make the link with previous investigations, we now 
consider the impact of a localized scatterer on the equipartition state. In Figure [3l we present a detailed 
balance of the energy carried by each wave mode to show that equipartition is globally preserved in a 
scattering system. Let us consider a bundle of plane P waves coming from space direction n' with mean 
squared amplitude ((\A P (n')| 2 )). Upon scattering, the P waves will be mode converted into S waves prop- 
agating into direction n. Reciprocally, we may consider incident S waves propagating in direction — n 
with mean-squared amplitude ((\A S (— n)| 2 )) and mode converted into P waves propagating in direction 
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Let us now consider the balance of energy carried by P waves in the scattering process: 
dE p 



dt 



d 2 n<(|^(n')| 2 }}(|/ 6l ^(n,n')| 2 + |/^(n, n')| 2 )dV 



d 2 n(((|^(-n)| 2 )}|r^(-n',-n)| 

7TX4-7T 

<<|^n)| 2 ))|/^(-n',-n)| 2 )dV, 



(94) 



where E p denotes the total energy carried by P waves. In equation (|94j) . the integrals are performed 
over all possible space directions n and n'. The first term represents the rate of energy conversion from 
P waves to S waves in the scattering process. This contributes to a decrease to the total P energy, hence 
the minus sign. The second term corresponds to the reciprocal process where P energy increases because 
of the S to P mode conversions. For the description of the shear wave motion, we have chosen a basis 
of linear polarizations {Si,^} as illustrated in Figure [2j Other bases such as left-handed and right- 
handed circular polarizations are popular in optics. For a wavefield at equipartition the mean squared 
amplitudes are independent of the vectors n and n' and satisfy the relations: c^{(\A Sl | 2 )) = Cp((|A p | 2 )), 
{(|j4 1 1 2 )) = {{\A S2 1 2 }}. Because of the fundamental reciprocity relations of conversion scattering derived 
in equation (|50|) , we find that the total energy carried by P waves is constant, i.e. dE p /dt = 0. This 
implies that equipartition is globally preserved by the introduction of a scatterer in the medium. This 
simple argument highlights again the link between Green's function reconstruction and equipartition. 
In this paper, the link between symmetry relations in scattering and the reconstructi on of Green's 



functi o n from random noise sourc es has been generaliz e d to t he case of elastic waves (see ISnieder et al 



( 2008 ): Snieder and Fleurvl ( 2010l ): Margerin and Sato ( 2011 ) for the acoustic case). Using the Dirac 
calculus, a number of previously known symmetry and reciprocity relations for the scattering amplitude 
have been derived in a systematic way ( VaratharajuluL Il977t iDassios et all . Il987t lAkl 1992). Based on 
the general Green's function reconstruction formula for stationary random noise sources, we have obtained 
a general form of optical theorem valid in both near and far-field, which incorporates the effect of a finite 
amount of absorption. We also show that the generalized optical theorem can be established from the 
scattering formalism in the limit of infinitesimal absorption. Thus, the correlation properties of random 
wavefields offer an alternative approach to study scattering phenomena in heterogeneous media. The 
reconstruction of Green's function in 3-D is illustrated with the simple case of a point-like heterogeneity 
embedded in an isotropic matrix. Based on a regularization of the return free space Green's function, we 
are able to sum the Born series for such a defect. Usual results of the formal boundary value problem for 
a small scatterer with an arbitrary density contrast are recovered in the low frequency limit. The point 
scatterer satisfies the generalized optical theorem, which in turn implies the reconstruction of Green's 
function in the presence of random distribution of noise sources. The point scattering model offers an 
attractive approach to tackle more complex multiple-scattering media. Such developments will be the 
topic of a separate publication. 
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